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Abstract 

We analyze geometrically non-linear isotropic elastic shells and prove the exis- 
tence of minimizers. In general, the model takes into account the effect of drilling 
rotations in shells. For the special case of shells without drilling rotations we present 
a representation theorem for the strain energy function. 

1 Introduction 

The paper is concerned with the geometrically non-linear 6-parameter resultant shell the- 
ory. This model of shells involves two independent kinematic fields: the translation vector 
field and the rotation tensor field, which have in total 6 independent scalar kinematic 
variables. This shell theory was originally proposed by Reissner [TT] and was developed 
consistently by several authors jHl |3] . 

In Section [2] we briefly present the kinematics of 6-parameter shells, as well as the equations 
of equilibrium and the constitutive equations of elastic shells. We formulate the boundary- 
value problem and prove the existence of minimizers associated to the deformation of 
isotropic shells in Section |3l This model is able to describe the effect of drilling rotations in 
shells. In order to see the difference to the Reissner-type kinematics for shells, we analyze 
in Section H] shells without drilling rotations and give a representation theorem for the 
strain energy function. Finally, we consider isotropic shells without drilling rotations and 
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identify the constitutive coefficients by comparison with the 6-parameter model. In case 
of shells without drilling rotations the strain energy function is only positive semi-definite, 
so that the general theorem for the existence of minimizers does not apply. 

2 Equations of equilibrium 

We consider a shell with the base surface in the reference cofiguration characterized by 
the position vector (relative to a fixed point O) : w C — )■ R^, = y^{xi,X2), and 
the structure tensor : w C — > 5*0(3), = d^{xi,X2) (X) ej . Here, (xi,X2) are the 
material curvilinear coordinates on , {ej} is the fixed orthonormal vector basis of the 
Euclidean space, and {d^{xi, X2)} is the orthonormal triad of directors which characterizes 
the orthogonal tensor field Q° |3i [5] . We employ the usual notations: the Latin indexes 
... take the values {1,2,3}, the Greek indexes ... the values {1,2}, the partial 
derivative daf = df /dxa , as well as the Einstein summation convention over repeated 
indexes. The set u is assumed to be a bounded open domain with Lipschitz boundary in 
the OxiX2 plane. 

For the deformed configuration of the shell, we denote by S the base surface, y{xi,X2) 
the position vector and {di{xi,X2)} the orthonormal triad of directors. The displacement 
vector is defined as usual hj u = y — y^ and the elastic rotation (between and S) by the 
proper orthogonal tensor field Q = di® . The orthogonal tensor R = QQ^ = di® Ci 
describes the total rotation from u to S. 

Let = day^ be the (covariant) base vectors in the tangent plane to 5°, n° = ai x 
02/11 tti X 02 II the unit normal to 5*°, and {«"} the reciprocal (contravariant) basis in the 
tangent plane, with ■ cip = and a" ■ = 5^ (the Kronecker symbol). Then the 
shell deformation gradient tensor is expressed by i*" = Grad^y = d^y ® a". 
We designate by N and M the internal surface stress resultant and stress couple tensors 
of the 1** Piola-Kirchhoff type for the shell, and by / and c the external surface resultant 
force and couple vectors applied to points of S, but measured per unit area of . The 
equilibrium equations for 6-parameter shells are [5] 

Div,Ar + / = 0, 

Div, M + ax\{NF^ - FN^) + c = 0, 

where Div^ is the surface divergence, {■)^ denotes the transpose, and axl( ■) represents the 
axial vector of a skew-symmetric tensor. We consider boundary conditions of the type 

m 

Nv = n*, Ml/ = m* along dSj , , 
y = y*, R = R* along dS^ , ^ ^ 

where u is the external unit normal vector to the boundary curve dS^ (lying in the tangent 
plane) and {dS'j , dS^ } is a disjoint partition of dS'^ . 
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According to [3] , the elastic shell strain tensor and the bending-curvature tensor 
K"^ in the material representation are 

= Q^F - Grad, = {Q^d^y - a„) ® a" 
= axl(g^9„Q) ® a". 



Under the hyperelasticity assumption, N and M are expressed by the constitutive equa- 
tions 

dW dW 

where W = W{E^, K^) is the strain energy density of the elastic shell. The boundary- 
value problem describing the deformation of non-linear elastic shells consists of the equa- 
tions We assume the existence of a function A{y,R) representing the potential 
of external surface loads /, c, and boundary loads n*, m* [4J. This loading potential can 
be decomposed additively as 

A{y, R) = Aso{y, R) + Aas^iv, R), 

Aso{y,R)= [ f-udS + Uso{R), 

^aso{y,R)= / n*-udl + U9so{R), 

where the load potential functions II^o , Ilg^o : /^^(w, 50(3)) M are assumed to be 
continuous and bounded operators. Corresponding to the deformation of elastic shells, 
we consider the following two-field minimization problem: find the pair (y, R) in the 
admissible set A which realizes the minimum of the functional 



Iiy,R)= WiE\K^)dS-Aiy,R) (5) 

J so 

for {y,R) e A, where A := {{y,R) G H\uj,R^) x H\co,SO{3)) \ y^Qs^ = y\ Rios^ = 
i?*}. Here, the boundary conditions are to be understood in the sense of traces, 
denotes as usual the Sobolev space, and L'^ the Lebesgue space. The variational principle 
associated to the total energy of elastic shells (jS]) has been proved in [3] . 

3 Shells with drilling rotations: existence of 
minimizers 

In case of physically linear isotropic shells, the strain energy density is assumed as the 
quadratic form 

2W{E'',K'') = ai{tTElY + a2tT{ElY + a3tT{E''^^'^El) +a4{n°Ey 
+/3i(trK|j)' + /32tr(Kf)' + /^g tr(K;'^K^) + (3^in'Ky, 
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where El = E^-n°^n^E^, Kjj = i<:'-n°®n°K'. The eight coefficients , Pk can 
depend in general on the initial structure curvature tensor = ax\(^daQ^Q^'^) ® 
For the sake of simplicity, we assume in our discussion that ak and f3k are constant. 

Theorem 1. Assume that the initial position vector is continuous and injective and 

do^y'^ e L°°{uj, R^), det (a„^(xi, xs)) > > 0, 

where qq 'is a constant. The external loads and boundary data are assumed to satisfy the 
conditions 

y*eH\iu,R^), R* e H\uj,S0{3)). 
If the constitutive coefficients satisfy the conditions 

2ai + 02 + as > 0, ^2 + > 0, 

- a2 > 0, a4 > 0, 2/3i + /32 + /^s > 0, (7) 
/32 + /33 > 0, /33 - /32 > 0, /34 > 0, 

then the minimization problem admits at least one minimizing solution pair {y, R) G 
A. 

Proof. In view of the inequalities ([7]) we can deduce that there exists a constants Cq > 
such that 

W{E'',K') > Co{\\E^f + \\K'f). 

Moreover, in this case the strain energy density W{E'^,K'^) is a strictly convex function 
of its arguments. Then, according to Theorem 6 from [1], we derive the existence of 
minimizers. The proof is based on the direct methods of the calculus of variations. □ 

Remark 2. For isotropic shells, the simplest expression of W{E^, K^) corresponds to 
the form ([6]) with 

«! = Cu, a2 = 0, as = C(l — u), 

a4 = asC(l -z/), Pi = Diy, /32 = 0, (8) 
/33 = D{1 - z/), /34 = atD{l - u), 

where h is the thickness of the shell, E the Young modulus, i> the Poisson ratio of the 
material, C = Eh/{1 — u'^) is the stretching (in-plane) stiffness, D = E h^/12{l — i/^) is 
the bending stiffness, and a^ , at are two shear correction factors. Note that the conditions 
([?]) are fulfilled for the coefficients ([8]). □ 

Without loss of generality, one can choose the directors {d^} such that = nP is the 
unit normal to S^. In what follows, we assume that ^3 = 71*^. 
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In the 6-parameter shell theory the drilling rotations are taken into account. The drilling 
rotation in a given point S can be interpreted as the rotation about the director ^3 . The 
general form of rotations about is 

Re = d3®d3 + cos6'(l-d3 (g) da) + sm9{dsxl), 

where 6 = 0{xi, X2) is the rotation angle and 1 = dj is the unit tensor. 

Let us describe next shells without drilling rotations in the framework of the 6-parameter 

shell theory. 



4 Shells without drilling rotations: 
characterization 

In case of shells without drilling rotations the strain energy density W must be insensible 
to the superposition of rotations Rg about d^ . This means that W{E'^ ^ K^) is assumed 
to remain invariant under the transformation 

Q ReQ. (9) 

In view of the definitions ([3]), this is equivalent to 

W{E', K') = W{[Q^R^day -aj® a", axl[Q^ RJ da{ReQ)] ® a") (10) 

for any angle 6{xi,X2)- The following result gives a characterization of shells without 
drilling rotation. 

Theorem 3. Assume that the strain energy function W is invariant under the trans- 
formation (Q . Then, W can be represented as a function of the arguments 

W = W{F^F , dsF , F^Grad.ds). (11) 

Conversely, any function W of the form f[TT]) is invariant under the superposition of 
drilling rotations ([9]). 

Proof. Firstly, it is clear that the function (ITT]) is invariant under the drill rotation, 
since F = Grad^y and Gradgds are both independent of rotations about d^ . Conversely, 
let us assume that W is invariant under the transformation (Q. If we denote by d^ = 
Redi = cos6'di + sin 6*^2 and dj — Red2 = —sinOdi + cos6d2 , then we find RgQ = 
d^ ® d^ . Inserting the last relation into equation f lTO|) and imposing the conditions that 
the derivative of ffTOj) with respect to 6 and 8^6 are zero, we obtain the equations 

dW dW 

_.e,..,„),_..,K.,K«).0 ^^^^ 
d{n^K ) 
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where we have used the notations a = aa a" = and c = d'j' ® cZg — ® . 
We interpret the relation flT^ i as a first order hnear partial differential equation for the 
unknown function W{E^,K'^), which depends on 12 independent scalar arguments (the 
12 components of E'^ and K"^ in the tensor basis {d^ ® cl"}). According to the theory 
of differential equations (see e.g., [I2], Chap. 6), to solve equation (|T2l) i we determine 11 
first integrals of the associated system of ordinary differential equations 

^ = c{E^ + a), ^ = c{K^ + K'). (13) 
ds ds 

We observe that the functions Uk are first integrals: 



Indeed, we have 



Ui = F^F = {E^ + af{E' + a), 
U2 = 'nPE" = dsF, Us = n^K\ 



(14) 



ds vds; .^^+«) + (^^'^+«)^-o, 

ATT J jpe 

^ = = n^\c{E' + a)] = 0, 

ds ds 

in view of f[T^ 1 , and analogously for C/3 . Finally, another independent first integral is 
the function 

Ui = F^Grad.dg = {E^ + afc (K^ + K^), (15) 

since ^ U4, = by virtue of relations f lT3|) . The functions f|T^ and f|T5|) represent in total 
11 scalar independent first integrals. Then, the general solution of the first order partial 
differential equation f lT^ i is 

W = W{F^F, dsF, n°K', F^Grad^dg), 

which in view of flT2|) 9 reduces to ffTTj) . □ 
Remark 4. From Theorem 3 follows that the strain energy f ITT]) can be alternatively 
expressed as a function of the following arguments 



* = (F^Grad.ds - Grad^n^) - S Grad^n^ ^ ^ 

= {E^^^^ a)cK^+ [ \ E^^^E^^ skew^;|j ]6, 

where we denote by b = — Grad^n'^. The tensor £ is a second order symmetric tensor 
accounting for extensional and in-plane shear strains, 7 is the vector of transverse shear 
deformation, and is a second order tensor for the bending and twist strains. □ 
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The results ( I16p are similar to those presented by Zhilin [T3j for shells without drilling 
rotations. The tensors S and 7 coincide with those given in [T3], but nevertheless the 
bending-twist tensor ^ is different. In this respect, Zhilin proposed the tensor 

= [F^ids X Grad.dg) + xb] + £{n^ xb). ^ ' 

We consider that the definition of the bending-twist tensor in the form f|T6|l ^ is more 
appropriate since the relation f|T7|) introduces an additional (unnecessary) rotation of 
Gradgds in the plane {£Zi,£Z2}- 

From f|T6|) and f|T7|) we see that in the linearized theory these deformation tensors reduce 
to: 

£ = sym(a Gradsw), 7 = n^Grad^u + ct/*, 
^ = c$ = cGiads{axp) + [skew(a Gradsu)]b, 

where if) is the vector of small rotations. One can easily see that £, 7 and 'J' are indepen- 
dent of the drilling rotation (if) ■ n^). In this case one gets the Reissner-type kinematics 
of shells [131 ISj with 5 degrees of freedom. 



5 Isotropic shells: comparison 

The isotropic shells without drilling rotations have been investigated in details by Zhilin 
[H], who determined the form of the strain anergy density ly as a quadratic function 
of its arguments (£,7,$). Suggested by these results, we consider the following strain 
energy function for elastic shells without drilling rotations (for the simplified case when 
the coefficients are independent of -K'") 

2Wi£n,^) = C[il-u)\\£r + i^itT£f] + |C(l-z/)/c7' no^ 
+D[i(l-z/)||*||2 + l(l-i/)tr(*2) + ^,(tr*)2], ^ 



where k is a shear correction factor. If we insert the expression f[T6l) of £,7 and ^ into 
f lTSjl . then we find the form of W in terms of the strain tensors (-E^, K'^). We observe that 
the resulting energy density W{E^, K^) is a super-quadratic function of its arguments. In 
the case of physically linear shells, when only the quadratic terms in {E^, K^) are kept, 
we obtain the simplified expression of the energy density 

+D[tT{Kl'^Kl) - ^ (trKf)2 -z/tr(Kf)2]. ^^^^ 



Finally, if we compare the relation f[T9|) with the general form of the strain energy density 
for isotropic shells ([6]) we find the following values for ak and 



ai = C z/, a2 = = C ^-y- , ai = C k. 



2 ' -4 - 2 (20) 
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Remark 5. The coefficients and given in f l2U]) for shells without drilling rotations 
are different from the values ([H]) corresponding to shells with drilling rotations. This 
indicates that the two types of shells will have different mechanical responses. □ 

Remark 6. The conditions ([7]) which insure the existence of minimizers are not satisfied 
by the coefficients (120]) corresponding to shells without drilling rotations since a^ — a2 = , 
2(3i + (^2 + = 0, and = 0. In this case, the strain anergy function (Q is not uniformly 
positive definite, and therefore the proof of existence of minimizers is more difficult (in 
this respect, see [8]). This is in accordance to the results presented by Neff [3 |8] for a 
plate model derived directly from the 3D equations of Cosserat elasticity. The comparison 
between the 6-parameter resultant shell theory and the model developed in [71 E] has been 
presented in [21 [1] . □ 
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